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Abstract 

We investigate the kinetic freeze out scenario of a nucleon gas through a finite 
layer. The in-medium mass modification of nucleons and it's impact on the freeze out 
' process is studied. A considerable modification of the thermodynamical parameters 

■ temperature, flow-velocity, energy density and particle density has been found in 

, comparison with evaluations which use a constant vacuum nucleon mass. 

\Q ■ 

1 Introduction 

As the fireball of a nucleus-nucleus collision cools down below a certain freeze out tem- 
perature Tpo < T c ~ 175 MeV the inelastic (chemical freeze out) and elastic (kinetic 
freeze out) collisions among the hadrons cease. This process is usually called the freeze 
out scenario (FO). Several approaches have been applied for the description of the freeze 
out of strongly interacting matter. Especially, there are kinetic models |2] as well as 
hydrodynamical approaches |3] which have been proven to be able to describe most of 
the collective phenomena like the different flow components in heavy-ion reactions. De- 
spite the success in comparison with experiments, the in-medium modifications of the 
hadrons during the freeze out process have not been taken into account yet. In all of 
the former evaluations the vacuum parameters of the particles have been implemented. 
However, during the freeze out process, the temperature and particle densities are presum- 
^ ably close to the deconfinement phase transition critical values jl]. Accordingly, strong 

in-medium modifications of hadronic properties like mass, width, coupling constants etc., 
are expected. The question arises how strong their impact on the freeze out process is. 
For answering this issue we study a nucleon gas and investigate how strong the impact 
of an in-medium nucleon mass shift on the thermodynamical parameters, temperature, 
flow-velocity, energy density and particle density, of the freeze out scenario is. 
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2 Freeze Out Process within a finite layer 

The theoretical description of the kinetic freeze out within a hydrodynamical approach 
has been worked out some years ago; see [5] and references therein. Here, the evaluations 
presented are based on the approach of Ref. where the freeze out description has 
been generalized to the case of a finite time-like layer (a finite space-like layer has also 
been investigated in more details about the approach presented will be published in 
|Sj). Accordingly, local equilibrium implies that the thermodynamical parameters inside 
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the layer become space-time dependent, i.e. we have a space-time dependent temperature 
T(x), flow velocity v(x), energy density e(x) and nucleon density n(x). From the equations 
of hydrodynamics we derive the following set of coupled differential equations for three 
unknowns 

de(x) = u tl {x)dT^ v {x)u v {x) + 2du lx {x)T' xv {x)u v , (1) 
dn(x) = u^{x)dN^(x) , (2) 

du»(x) = -r-r (g^ - u^xWtx)) dN v . (3) 
n(x) \ J 

The differentials in © - © are deduced from the microscopic definition, dN^(x) = 
J Cr& M df(x, k) for the particle current, and dT fMU (x) = J ^k^k u df(x, k) for the energy 
momentum tensor. Here, x M = (t, r) is the four-coordinate and fc M = [E^,\s) is the four- 
momentum of the nucleon; the scalar one-particle distribution function is normalized to 
the invariant number of nucleons N, i.e. N = J d 3 r c/ 3 k f(x,k). We mention that the 
second term in © vanishes within the approach presented, but not in general. The 
Rest Frame of the Gas (RFG) is defined by u^ FG (x) = (1, 0, 0, 0). For baryon-dominated 
matter, like in our case, the Eckart's definition jHj of four- flow is commonly used, defined 
by u^i^x) = N^(x) I ^N v {x)N v {x). Another Lorentz frame, e.g. Rest Frame of the Front 
(RFF), can be defined by a Lorentz boost in respect to RFG. 

Since there are four unknowns in the problem under consideration an additional con- 
straint is necessary, which is provided by the Equation of State (EOS) for the nucleon gas 
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e(x) = n(x) 



M N (n(x),T(x)) -E + ^ - 1 ) + jj T{x) 
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(4) 



The nuclear binding energy is Eq = 16 MeV, and K =~ 235 MeV is the compressibility; 
M N (n,T) is the in-medium nucleon pole mass. The EOS (HJ) is used to determine the 
temperature T{x) of the interacting component of the nucleon gas during the freeze out 
process. Accordingly, the four equations ©, (J2J), © an d © represent a closed set for 
evaluating the four unknowns T, v, e, n of the one-particle system. 

Furthermore, as the system expands and cools down the number of interacting par- 
ticles decreases up to the post freeze out surface of the finite layer, where by definition 
the density of interacting particles vanishes. Correspondingly, the particle distribution 
function is decomposed into two components, an interacting part f\ and a non-interacting 
part ff, thus f(x,k) = fi(x,k) + ff(x,k). Accordingly, there is an interacting particle 
density rij and a non- interacting particle density rif with n = rif + On the pre- freeze 
out hypersurface we assume to have thermal equilibrium, i.e. we have a Jiittner distri- 
bution, cf. 0, for fi as starting one-particle distribution function, while by definition ff 
is zero on the pre-freeze out hyper-surface. The space-time evolution of the interacting 
and non-interacting components inside the layer is governed by the following differential 
equations 

*'< - -;(i^)(^)^ + > (t) - /J " (5) 
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with the time r between collisions, and f eq is the Juttner distribution, cf. ||. The second 
term in (jSJ) is the re-thermalization term [S], which describes how fast the interacting 
component approaches the Juttner distribution within a relaxation time Tq. Here, we will 
use the immediate re-thermalization limit tq — > 0, which implies fi —>■ f eq faster than 
ro — > 0. The explicit expressions for the differentials dN^ and dT^ u within the approach 
presented can be found in [Hj. 

The set of equations (0) - © allow us to evaluate the basic thermodynamical function 
T(x), v (x), e(x) and n(x) during the freeze out process for a particle with a mass M N (n, T). 



3 Nucleon Mass Shift 

A generalization of the QCD sum rule approach for nucleons at finite densities and tem- 
peratures ^T] leads to the following expression for the pole mass of a nucleon embedded 
in a hot and dense hadronic medium, 

M N (n,T) = M N (0)+ReX s (n,T) + ReX v (n,T) , (7) 

where M^(0) = 939 MeV is the vacuum nucleon pole mass, and with the attractive scalar 
part (Re S 5 < 0) and the repulsive vector part (Re Sy > 0) of nucleon self energy in 
medium. Typical values at nuclear saturation density = 0.17fm -3 and at vanishing 
temperature are Re S 5 = -400 MeV, Re £y = +300 MeV |TT|. Here, we will take the 
QCD sum rule results for a nucleon in matter, given by ^1] 

ReX s (n,T) = +M K (0) - l) , (8) 




Re £y(n, T) = ~^<0)™ • (9) 

Here, (Q\qq\Q) and (Q\q^q\Q) are in-medium condensates and \ is a state which describes 
the hot and dense hadronic matter inside the layer, while (0|gg|0) = (—0.250 GeV) 3 is 
the chiral condensate. In (JHJ) and © we have neglected the gluon condensate and higher 
mass dimension condensates which give rise to small corrections only. 

There are two nucleonic components inside the finite layer: an interacting component 
with density n« and a non-interacting component with density nf. For evaluating the 
condensates (jHJ) and © we approximate the interacting component by a Fermi gas with 
chemical potential /ij and temperature T. On the other side, the temperature for the non- 
interacting component becomes ill-defined. Nonetheless, a relevant physical parameter 
for describing the non-interacting component remains the density nf. Accordingly, the 
condensates in one-particle approximation are given as follows [TB"] : 

<n|*?|n> = (0|*g|0) + i(^,T) (N(k)\qq\N(k)) + -^l^(N(k)\qq\N(k)) , (10) 

(n|fftg|n> = j(^,t) (iv(k)|gtg|iv(k)) +^5j(^( k )l^l^( k )>> (ii) 

with I = 4 / d 3 k/[(2Ek) (27r) 3 ) (exp((Ek — fJ-ij/T) + 1)]; for vanishing temperature we have 
7 rii/(2M N ). Note that <0|g + g|0> = 0; the nucleon energy is E k = J M N (0) 2 + k 2 . For 
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nucleons we take the relativistic normalization (iV(ki)|iV(k2)) = 2E/ < . 1 (2ti) 3 S^ (ki — k 2 ) 
is used. The chemical potential for the interacting component can be evaluated via rii = 
4 / d 3 h/[(2n) 3 (exp((Ek — /ij)/T) + 1)]. The condensates in Fermi gas approximation are 
given by [12] (iV(k)|gg|iV(k)) = M N (0)a N /m q and (N(k)\q^q\N(k)} = 3M N (0). The 
nucleon sigma term is o~n — 50 MeV, and m q ~ 5 MeV is the averaged current quark mass 
of the light quarks. Inserting these parameters into (JHJ) and (jHJ) we obtain Re = —390 
MeV and Re £y = 315 MeV at ground state saturation density n . The Eqn. flj|) - 
(fTTj) summarize our propositions made for obtaining the in-medium nucleon pole mass 
M N (n,T), which enters the EOS ® and the differentials dN>* and dT^. 

4 Results and Discussion 

For all of the calculations, we have taken T preFO = 150 MeV, n pX ePO — 1-5 n and 
Vpve fo = 0.5 c as starting values on the pre freeze out hypersurface. These values 
are, for instance, in line with typical parameters which have been reached within the 
Alternating-Gradient Synchrotron (AGS) at Brookhaven National Laboratory (BNL) in 
Brookhaven/USA. Higher baryonic densities can be reached within the Schwer-Ionen- 
Synchrotron (SIS) at Gesellschaft fur Schwerionenforschung (GSI) in Darmstadt/Germany. 
Note that T pre fo and n prc fo are pre freeze out values and, therefore, they are larger than 
typical post freeze out values given, for instance, in Ref. jl]. 

In Figs. [Dand|2l the time evolution of the primary thermodynamical functions through 
the finite freeze out layer are shown, in terms of the proper time r. Note that the densities 
n = rii + rif and e = + Cf are kept constant inside the layer. 




Figure 1: LEFT: The temperature of the interacting component. RIGHT: The flow 
velocity parameter v of the interacting component. The solid lines are with a constant 
nucleon mass M^(0) = 939 MeV, while the dashed curves are evaluated with a density 
and temperature dependent nucleon mass Mn(u,T). 
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Figure 2: LEFT: Nucleon energy density of the interacting component. RIGHT: Nu- 
cleon particle density of the interacting component. The solid lines are with a constant 
nucleon mass Mjv(O) = 939 MeV, while the dashed curves are evaluated with a density 
and temperature dependent nucleon mass M N (n,T). 

We find a substantial impact of in-medium mass modification on the freeze out process 
within the purely nucleon gas model. Furthermore, the Figs. and El also elucidate, 
that the freeze out process proceeeds faster for all thermodynamical quantities T, v, e, n 
when taking into account the mass dropping of nucleons. The physical reason for a faster 
freeze out originates from a smaller energy density of the nucleon system due to a smaller 
nucleon mass M^(n,T) compared to the vacuum nucleon mass Mjv(O). 

Finally, we remark that in-medium modifications have actually to be taken into ac- 
count not only during the freeze out process, but also before, i.e. during the hadronization. 
This points then to an even stronger impact of in-medium modifications on the final par- 
ticle spectrum than presented. 



5 Summary 

We have investigated a freeze out scenario within a finite layer for a massive nucleon gas. 
Special attention has been drawn about the issue how strong the impact of the in-medium 
nucleon mass modification on the thermal freeze out process is. By focussing on a purely 
nucleon gas we have found a substantial effect on the thermodynamical quantities like 
temperature T, flow velocity v, particle density n and energy density e of the interacting 
component. All of these thermodynamical functions have revealed a faster freeze out 
compared to a scenario without an n-medium nucleon mass shift. 

In summary, our findings for a nucleon gas suggest that taking into account in-medium 
modifications of nucleons seems to be necessary and an interesting phenomenon, in par- 
ticular for collision scenarios with high baryonic densities. 
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